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Evidence for Astrophysical Muon Neutrinos from the Northern Sky with IceCube 
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Results from the IceCube Neutrino Observatory have recently provided compelling evidence for 
the existence of a high energy astrophysical neutrino flux utilizing a dominantly Southern Hemi¬ 
sphere dataset consisting primarily of and Vt charged current and neutral current (cascade) 
neutrino interactions. In the analysis presented here, a data sample of approximately 35,000 
muon neutrinos from the Northern sky was extracted from data taken during 659.5 days of live¬ 
time recorded between May 2010 and May 2012. While this sample is composed primarily of 
neutrinos produced by cosmic ray interactions in the Earth’s atmosphere, the highest energy 
events are inconsistent with a hypothesis of solely terrestrial origin at 3.7(7 significance. These 
neutrinos can, however, be explained by an astrophysical flux per neutrino flavor at a level of 
^{E,j) = O.Olg ® X 10“^® GeV“^ cm“^ sr“^ s“^ (consistent with IceCube’s Southern 
Hemisphere dominated result. Additionally, a fit for an astrophysical flux with an arbitrary spectral 
index was performed. We find a spectral index of 2.2^0 2 , which is also in good agreement with the 
Southern Hemisphere result. 


The nature of the objects and the mechanisms which accelerate cosmic rays pose major open questions in cur- 



3 


rent astrophysics, which may, in part, be answered by 
observations of high energy neutrinos. At high energies, 
the majority of cosmic rays are protons or atomic nuclei, 
and their interaction with other matter or radiation is 
known to produce neutrinos [T] . If this happens near the 
source of the cosmic rays, the neutrinos, which—unlike 
the charged cosmic rays—can travel undeflected through 
the magnetic fields of deep space, can point back to these 
sources. 

IceCube is a detector constructed at depths between 
1.5 km and 2.5 km in glacial ice at the South Pole, in¬ 
strumenting about a cubic kilometer of volume with op¬ 
tical sensors [2] • This forms a Cherenkov detector for the 
light produced when neutrinos interact and generate sec¬ 
ondary charged particles. These interactions give rise to 
two characteristic event topologies: linear ‘tracks,’ pro¬ 
duced by long-range muons emitting light as they travel, 
and near-spherical ‘cascades,’ from the more point-like 
light emission of electromagnetic and hadronic particle 
showers which terminate in ice after small distances com¬ 
pared to the instrumentation density of the detector [3] . 

One effective method for identifying neutrino interac¬ 
tions is to look for events which show no sign of light 
emission when entering the detector boundary. These 
are referred to as ‘starting’ events. A recent IceCube 
study using this technique [i] has determined that as- 
trophysical neutrinos at high energies do exist, and that 
their flux is broadly compatible with existing models 0 - 
[7]. While such starting events provide good evidence for 
an astrophysical neutrino flux, they do not sample all 
components of the expected flux equally well. Due to 
absorption in the Earth, few neutrinos are observed from 
the Northern sky, and few of the observed events are 
identifiably This analysis seeks to observe more of 
these particular types of events by relaxing the require¬ 
ment that events begin inside the detector to permit the 
use of the long muon range to achieve a larger effective 
volume. Events are then selected based on the event 
topology of muons produced from interactions to re¬ 
duce background contamination. In this analysis, as in 
other IceCube analyses, it is not possible to distinguish 
neutrinos from antineutrinos, so only the combined flux 
can be measured. 

To identify astrophysical muon neutrinos, we must dis¬ 
tinguish them both from other types of events in the 
detector and from other sources of neutrinos. The ma¬ 
jority of the data recorded by IceCube are produced by 
muons originating in cosmic ray air showers that pene¬ 
trate the ice and reach the detector. Since this analysis 
seeks to take advantage of the long muon tracks and can¬ 
not depend on observing the neutrino interaction vertex 
inside the detector, only muons with directions that im¬ 
ply they passed through more material than the maxi¬ 
mal expected muon range are selected. In this case, part 
of the distance must have been traversed by a neutrino, 
which is less prone to interaction. This analysis accepts 


therefore only events whose reconstructed zenith angles 
are greater than 85°, corresponding to an overburden 
equivalent to at least 12 km of water. The directions 
of muon events are reconstructed by fitting the hypoth¬ 
esis of a particle moving at the speed of light and emit¬ 
ting Cherenkov radiation to the timing of the observed 
photons. The fit accounts for the expected delay of the 
first photon to reach each detector module due to scat¬ 
tering [B]. Rejecting poorly fit events removes both low 
energy atmospheric muons with poor direction resolution 
and the less numerous cascade-like events produced by 
neutrino interactions other than charged-current In 
addition to the direction of the muon, the other observ¬ 
able of interest is muon energy. A proxy for the energy 
is computed by fitting the amount of light expected to 
be emitted by a template muon to the number of ob¬ 
served photons in each event m [ 10 ]. The precision of 
the energy proxy is limited by the relatively short sec¬ 
tion of the muon’s total track which is observed, and is 
only loosely connected to the energy of the interacting 
neutrino since an unknown amount of energy is gener¬ 
ally lost before the muon reaches the detector. After 
applying event-quality criteria (which are qualitatively 
equivalent to those used in earlier studies nnini, with 
details being given in the online supplement [13] and in 
[Tl] f this yields a highly pure (99.9%) sample of neutrino- 
induced muon events, with an efficiency of about 24% for 
neutrino-induced events from an E~^ spectrum. This se¬ 
lection still suffers from neutrino absorption in the Earth, 
resulting in a loss of events at the highest zenith angles 
and energies. This analysis was performed with a blind¬ 
ness criterion such that only 10 % of the experimental 
data were used in its development, in conjunction with 
simulated data, to determine the data selection The full 
data are used only after the analysis technique had been 
fixed. 

Since the astrophysical neutrinos we seek to observe 
in this study are expected to be produced in conjunc¬ 
tion with the cosmic rays [mill], they should have a 
related power-law spectrum of the form $ oc E~'^, where 
7 should be ^ 2. For this analysis we take 7 = 2 as a 
benchmark model We also make the further simpli¬ 
fying assumption that the astrophysical flux is isotropic, 
as would be the case for a signal originating from many 
distant, individually weak sources. 

Although astrophysical neutrinos are the target of the 
analysis, the numerous neutrinos produced by cosmic ray 
air showers must be accounted for. Atmospheric neutri¬ 
nos are usually separated into two groups: those pro¬ 
duced by the decays of pions and kaons, referred to as 
‘conventional,’ and those produced by the decays of heav¬ 
ier mesons, particularly those containing charm quarks, 
referred to as ‘prompt’. Since the conventional atmo¬ 
spheric neutrinos arise from relatively well-understood 
particle physics and have been measured by a variety of 
experiments [IHIIISI, there exist several models for this 
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flux Here we use the HKKMS07 calculation [5D], 

where the uncertainty of this calculation is estimated by 
its authors to be less than 10% at few GeV energies, 
which is consistent with measurements [23) . and is ex¬ 
pected to increase with energy to around 25% at 1 TeV. 
Since this model was designed for relatively low ener¬ 
gies (100 MeV-10 TeV) compared to those considered in 
this analysis (^100 GeV-100 TeV), it is extended and 
modified according to the procedure in m to take into 
account the input cosmic ray spectrum |24j at high en¬ 
ergies. An important feature of the conventional atmo¬ 
spheric neutrino flux is that the parent mesons may be 
destroyed by interactions with the medium before decay¬ 
ing and producing neutrinos. The energy spectrum is 
therefore steeper (oc than that of the cosmic rays 

from which it is produced (oc [^. This is then 

markedly softer than the hypothesized spectrum of as- 
trophysical neutrinos. The cosmic ray showering process 
gives these neutrinos a characteristic distribution in di¬ 
rection, peaked near the observer’s horizon, because of 
the different profiles of atmospheric density the air show¬ 
ers encounter. 

The prompt atmospheric neutrinos are less well un¬ 
derstood, as they have not yet been observed experimen¬ 
tally, and the theoretical predictions depend on under¬ 
standing heavy quark production in cosmic ray-air col¬ 
lisions at high energies. Multiple calculations exist Hi- 
[2H]; and here we choose the phenomenological ERS esti¬ 
mate of the flux [3S], again applying corrections for the 
input cosmic ray spectrum. This model has a normal¬ 
ization uncertainty of about a factor of two, and other 
calculations predict substantially larger or smaller fluxes. 
Like the conventional atmospheric neutrinos, the energy 
spectrum of the prompt component arises from the spec¬ 
trum of the cosmic rays. However, since the intermediate 
mesons involved decay so rapidly (with a mean lifetime of 
1.04x 10“^^ s for the at rest, as opposed to 2.60x 10“® 
s for the tt^ or 1.24 x 10“® s for the K^), losses via inter¬ 
actions are suppressed and the spectrum remains similar 
to E~'^-'^, and likewise remains essentially isotropic. 

To fit the observed data, we implement the binned 
Poisson profile likelihood construction described in m- 
Here, the expected event rates for each flux component 
are computed by weighting a generalized simulation of 
neutrinos traversing the Earth and interacting at IceGube 
according to the model’s input neutrino flux. Gompar- 
isons are made in each bin to the observed data. For 
this study, the data are binned in both the reconstructed 
zenith angle and the energy proxy. The main parameter 
of interest for this fit is the normalization assigned to the 
astrophysical flux component, while the normalizations 
of the background components are treated as nuisance 
parameters. Additional nuisance parameters include the 
difference between the true slope of the cosmic ray spec¬ 
trum and the assumed model, the efficiency with which 
the IceGube hardware detects photons emitted in the ice. 



FIG. 1. The distribution of reconstructed zenith angles of 
events in the final sample, compared to the expected distribu¬ 
tions for the fit of an E~^ astrophysical neutrino spectrum. 
Only statistical errors are shown, though in almost all bins 
they are small enough to be hidden by the data markers. 



FIG. 2. The distribution of reconstructed muon energy proxy 
for events in the final sample, compared to the expected distri¬ 
butions for the fit of an E~’^ astrophysical neutrino spectrum. 
Only statistical errors are shown. The energy proxy does not 
have a linear relationship to actual muon energy, but values 
~ 3 x 10^ are roughly equivalent to the same quantity in GeV. 
Larger proxy values increasingly tend to underestimate muon 
energies, while smaller values tend to overestimate. 

and the relative contributions to the conventional atmo¬ 
spheric neutrino flux from kaon decays rather than pion 
decays. The nuisance parameters can be constrained us¬ 
ing prior information from external sources, and the pri¬ 
ors used in this analysis are fisted in the fourth column 
of Table HI 

The parameter values from fitting 659.5 days of de¬ 
tector livetime using the benchmark set of fluxes are 
summarized in Tab. |lj and the projections of the ob¬ 
served and fitted spectra into the reconstructed zenith 
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Parameter 

E-^ Fit 

Best Fit 

Prior 

Astrophysical flux normalization per flavor Q.Olg ® x 10 
Astrophysical flux index fixed to 2 

1.7tg j X 10"^® 

c\ rt+0.2 
^•^-0.2 

> 0 

none 

HKKMS07 normalization 

ERS normalization 

Cosmic ray spectral index change 
Detector optical efficiency 

Kaon production normalization 

o.93l°;o! 

-0.024«°“ 

+9.ll°;®% 

1 1 r+O.OS 

07 

n qq+0 04 
'^•^'^-0.04 
q+1.05 

-0.023tZs 

+9.lt°oi% 

1 1 i^+0-08 

> 0 

> 0 

0T0.05 

none 

liO.l 


TABLE I. Fit parameters are shown for two case: when an E~^ astrophysical flux with equal flavor composition and equal 
neutrino and antineutrino components is assumed {E~^ Fit), and when the index of the astrophysical flux is allowed to vary 
(Best Fit). The listed error ranges are 68% confidence intervals. The gaussian priors are shown as the mean value ± the 
standard deviation, but the fit results do not depend substantially on the priors. Units for the astrophysical flux normalization 
are GeV”'^ cm“^ sr“^ s“^, and HKKMS07 |20| and ERS are the reference conventional and prompt atmospheric fluxes, 
respectively. 


angle and muon energy proxy are shown in Fig. [l] and 
Fig. H respectively. The uncertainties shown for the fit 
parameters include both statistical and systematic con¬ 
tributions (at the 68% confidence level), via the profile 
likelihood, using the approximation [2H]. Note that 
the data point in Fig. [^at muon energy proxy values of 
around 1.4 x 10^ should not be taken as an indication of a 
spectral feature: A fluctuation of this size is expected to 
occur in approximately 9% of experiments due to statis¬ 
tical fluctuations, and even a delta function component 
in the true neutrino spectrum would be broadened into 
a far wider peak in the muon energy proxy [10) . 

The best fit for the astrophysical com¬ 
ponent is a flux $(i?j,) = 9.oil;® x 

10~^^ GeV~^ cm~^ sr~^ s~^ ( ^gi^ev ) ^ flavor. 

The best fit prompt component is 0.94 times the bench¬ 
mark flux, but is consistent with zero. The significance 
of the non-zero astrophysical flux is evaluated by a 
likelihood ratio test to the null hypothesis that only 
atmospheric neutrino fluxes are present, in which case 
the fitted prompt atmospheric normalization rises to 
4.0 times the ERS model. An ensemble of trials is used 
to establish the distribution of the likelihood ratio test 
statistic, yielding a p-value of 1.1 x 10“"^ or a single-sided 
significance of 3.7cr. 

The range of neutrino energies in which this astrophys¬ 
ical flux is constrained by the data is calculated to be 330 
TeV-1.4 PeV. The endpoints of this range are found by 
applying a hard cutoff to one end of the astrophysical flux 
template, refitting the data with the other astrophysical 
flux parameters held constant, moving the cutoff inward 
until the resulting fit likelihood is 0.5cr worse than the 
best fit. This gives a conservative estimate of the energy 
range in which the astrophysical flux is necessary to ex¬ 
plain the observed data, although the flux may actually 
have a greater extent [30]. The flux should not be inter¬ 
preted as existing strictly within this energy range; were 
this the case simulation trials suggest that this analysis 
would measure a flux normalization only 5-20% of the 



Power law index 


FIG. 3. Likelihood profile of the astrophysical flux power- 
law index and the flux normalization at 100 TeV in units 
of 10“^® GeV“^ cm“^ sr“^ s“^. While the E~‘^ result is well 
within the 68% contour, it is not the overall best fit. Also 
shown are the best fits from various IceCube analyses of start¬ 
ing events, which generally have good agreement: Starting 
Events (HE) |4|, Starting Events (LE 1) [31], Starting Events 
(LE 2) [31]. 


result shown in Table H] 

Since the true flux need not have a spectral index 
of exactly 2, the fit was repeated allowing the index 
to vary, leading to a result of ^{Ey) = 1.7lQ;g x 

10~^^ GeV~^ cm~^ sr~^ s~^ ( lo^ev ) The nui¬ 

sance parameters do not change significantly except the 
prompt atmospheric normalization, which falls to zero, as 
shown in Tab. [I] Figure[^shows the confidence regions for 
the astrophysical flux normalization and spectral index, 
and compares this result to three other IceCube analy¬ 
ses using starting events [HEIlEi]. The compatibility 
of these results is noteworthy because this work uses an 
independent set of data from the others (a single, near¬ 
horizontal, high energy track event is shared with the 
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other samples), while the starting event results are highly 
correlated with each other. The spectral indices found by 
this work and by the starting event analyses are consis¬ 
tent within their respective uncertainties, but the best fit 
spectrum for this data set is slightly harder than those 
for the starting event analyses, particularly those extend¬ 
ing to lower energies, which are uniquely able to probe 
non-atmospheric contributions to the neutrino flux. A 
single power law provides an acceptable fit to all data, 
however, the present data cannot yet rule out the possi¬ 
bility that the astrophysical neutrino flux is not isotropic 
or that the spectrum is not a pure power law. 

In this study we see a clear excess of data above the 
expected atmospheric neutrino backgrounds at high ener¬ 
gies, similar to the result of [4j. In particular, despite the 
fact that these are almost entirely disjoint datasets (a sin¬ 
gle, near-horizontal track event, event 5 from [J], appears 
in both samples), both excesses are consistent in nor¬ 
malization within uncertainties, assuming an E~^ spec¬ 
trum: 9.5±3x 10“^®GeV“^ cm“^ sr“^ s“^ from the start¬ 
ing event study and 9 . 9^3 4 x 10 “^®GeV“^ cm“^ sr“^ s“^ 
from this work. These measurements do use different 
calculations of the neutrino-nucleon cross-sections, which 
influence the conversion of the flux into a rate of observed 
events: The starting event study used the calculation of 
[55] . while this study uses the updated calculation from 
[55] . which differs by 5-10% at the energies relevant to 
these analyses, but this is a relatively small effect com¬ 
pared to the uncertainties of these results. Thus, the 
observed data are found to be consistent with a flux 
consisting of equal parts of all neutrino flavors. Simi¬ 
lar consistency is seen in a recent analysis of starting 
events [52] • As shown in Fig. the results for arbitrary 
power laws are also in good agreement. These two mea¬ 
surements are compared in Fig. [^ along with other re¬ 
cent measurements and theoretical models. The result of 
this study also suggests that astrophysical neutrinos are 
present at the several hundred TeV energies where ob¬ 
servations were lacking in the dataset of | 3 ], suggesting 
that this was merely a statistical fluctuation. 

Models of the astrophysical neutrino flux besides un¬ 
broken power laws can also be considered. Here we ex¬ 
amine a small number of representative models. One 
candidate source type is the cores of active galactic nu¬ 
clei (AGN) [51 I55H55] . A fit of the AGN flux model |5] to 
the data in this analysis demonstrates in an incompati¬ 
bility in the normalization, with the predicted flux being 
too large by a factor of 6 . Another possible source class 
are regions with high star formation including Starburst 
galaxies [51 I55H55] . Gomparing the spectrum pro¬ 

posed by [5] to the data reported here, we find that it 
is compatible after its normalization is multiplied by a 
factor of 2.5. Finally gamma ray bursts (GRBs) have 
been long considered candidates for neutrino production 
[ZmiHlT], but recent dedicated searches by IceCube for 
neutrinos correlated with GRBs have placed strong limits 



FIG. 4. Comparison of the best fit per-flavor astrophysical 
flux spectrum of E~^ from this work, assuming a flavor ratio 
of 1:1:1, (shown in dark green with the 68% error range in 
lighter green) to other selected IceCube measurements (heavy 
lines) [U [ 12 ] and theoretical model predictions (thin, dashed 
lines) jlHZl 11711201128j . The sensitivity of this analysis is also 
shown as the thin, green line. 

disfavoring this hypothesis [55] . 

While this work represents the first strong evidence for 
an astrophysical Vfj, flux in the Northern Hemisphere, the 
sources producing these neutrinos remain unknown. Al¬ 
though muon events in IceCube have sub-degree angular 
resolution, recent IceCube searches for point-like and ex¬ 
tended sources of muon neutrinos found no statistically- 
significant evidence for event clustering, nor correlation 
of neutrinos with known astrophysical objects [59] . In the 
Northern Hemisphere, the point source flux upper limits 
are 10 — 100 times lower than the total diffuse flux level 
observed here, so the flux cannot originate from a small 
number of sources without violating those limits. The 
constraint on the number of sources was explored with a 
simple simulation where sources were injected uniformly 
over the Northern sky, with fluxes at the maximum levels 
allowed by the point source upper limit at each selected 
point, until the total flux reached the measured diffuse 
flux. On average, at least 70 sources are required to main¬ 
tain consistency with the point source upper limits. This 
assumes each source is a true point source and emits an 
unbroken E~^ ‘^ power-law flux. If the sources instead 
follow harder E~‘^ power law spectra, the diffuse flux 
could be split across an average of ~ 40 sources while re¬ 
maining consistent with the point source analysis. Given 
that the diffuse flux in the Southern Hemisphere is ob¬ 
served at a similar flux level, this observation suggests 
that the flux has a large isotropic component dominated 
by a large population of extragalactic sources, although 
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local sources can still have significant contributions. 
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DATA SELECTION 

This section describes in greater detail the quality cri¬ 
teria applied to select the events used in this study. These 
criteria were designed by comparing 805.5 hours of exper¬ 
imental data from the 2010 data-taking period (repre¬ 
senting 10% of the data from that period and 5% of the 
total data used in the full analysis) to simulated data. 
The simulation included a cosmic-ray dataset equivalent 
to approximately 264 hours of air shower background 
and a neutrino dataset weighted to both a conventional 
atmospheric neutrino spectrum and a hypothetical E~'^ 
powerlaw spectrum (representing a possible astrophysical 
flux). 

1. Muon Filter: The data available as input for this 
analysis are those selected by the IceCube detec¬ 
tor’s ‘online’ filters, which determine the subset of 
the data to be transmitted to the Northern Hemi¬ 
sphere. The filters operate on the output of the 
detector data acquisition system, which generates 
event readouts whenever a trigger condition is met. 
Several triggers are implemented, but a represen¬ 
tative example is the simple multiplicity trigger, 
whose condition is that some minimum number of 
detector modules (typically 8) observe light within 
a short time window (typically 5 microseconds), 
where each of the contributing modules must have 
made its observation within one microsecond of one 
of its four vertically neighboring modules. The trig¬ 
gers are designed to read out the detector when use¬ 
ful physics signals are likely to be present, which the 
filters then attempt to reconstruct and classify. The 
‘Muon Filter’ is designed to select events compat¬ 
ible with neutrino-induced muons passing through 
the detector. This includes upwards-going events 
as well as high-energy downwards-going events. 
Low-energy downwards-going events are mostly air 
shower muons which are uninteresting for this anal¬ 
ysis. Event energies are estimated very simplisti- 
cally at this level, using the total amount of ob¬ 
served light as a proxy. 


2. Data Reduction: Detailed reconstructions are per¬ 
formed on the data selected by the Muon Filter, 
which can reveal that the event should not have 
originally passed the filter criteria, so reapplication 
of the filter criteria can remove many of the air 
shower muons which are uninteresting to this anal¬ 
ysis. Additionally, a correct track reconstruction 
should generally pass closer to the points where 
more photons were detected, so the average dis¬ 
tance from the reconstructed track to the detec¬ 
tion points (the detector modules which observed 
light), weighted by number of photons detected, 
should be small. Events with up-going reconstruc¬ 
tions whose weighted distance was greater than 200 
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meters were therefore rejected, unless they were 
reasonably bright (more than 100 photoelectrons 
observed). Some of the air shower events are actu¬ 
ally coincident, with bundles of muons from unre¬ 
lated showers passing through the detector during 
the same readout window. A separation algorithm 
is also applied to attempt to separate these into 
distinct sub-events, which are then reconstructed 
normally. 

3. Zenith Angle: The main power of this selection to 
reject air shower muon bundles derives from select¬ 
ing events which pass through considerable over¬ 
burden to reach the detector. At some zenith an¬ 
gle, 0, the overburden is no longer sufficient and 
air shower events will reach the detector at a non- 
negligible rate. Events with cos(0) <0.1 (0 « 85°) 
must pass through at least 12 km of water equiv¬ 
alent to reach the top of the IceCube detector, so 
only these events are selected for use in this analy¬ 
sis. 

4. Track Quality: In most of the zenith angle range 
used by this analysis, the Earth blocks cosmic ray 
muons entirely from reaching the detector. Some 
such events are present in this part of the observ¬ 
able space in the reconstructed data, but only be¬ 
cause their reconstructions are inaccurate. Like¬ 
wise, neutrino-induced events whose directions are 
severely misreconstructed need to be rejected. Two 
methods are used here, both of which derive from 
the fact that the fit used to reconstruct these data 
is a maximum likelihood calculation, and addi¬ 
tional information can be extracted from the like¬ 
lihood description besides the parameter values of 
the best fit. The first variable used, known as the 
‘paraboloid sigma’ uses the shape of the profile like¬ 
lihood around the best-fit point to estimate the sta¬ 
tistical uncertainty on the location of that point; 
the second, the ‘reduced log likelihood’ (RLogL) 
attempts to use the best obtained likelihood value 
as a global measure of the success of the fit. The 
combined cut on these variables is shown in Suppl. 
Fig. [2 

5. Bayesian Reconstruction: Another means of elim¬ 
inating misreconstructed events is to compare the 
unconstrained track reconstruction with one which 
has incorporated a bayesian prior that the majority 
of observed events are truly down-going, and should 
be reconstructed as such. In the case that an event 
is correctly reconstructed, the evidence due to the 
data should overwhelm this prior, and both recon¬ 
structions will have similar results. In the up-going 
region, one can make the simplifying assumption 
that if the unconstrained reconstruction performed 
better by some constant factor (corresponding to 
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SUPPL. FIG. 1. The distributions of simulated air shower 
background events and E~^ signal neutrino events in two 
measures of directional reconstruction quality, whose mean¬ 
ings are described in the text (Data Selection item 4l. Many 
air shower events remain in the sample at this stage only 
because their reconstructed directions are incorrect, so these 
variables allow separating these types of data. 


the floor in the prior) than the constrained recon¬ 
struction, it was probably trustworthy. For events 
in the down-going region which are correctly recon¬ 
structed, the reconstructions with and without the 
prior should be similar, and thus the difference in 
their log likelihoods will be related by the prior, so 
the distribution of well-reconstructed events should 
approximately trace out the prior as a function of 
zenith angle. The cut which was used on this vari¬ 
able is shown in Suppl. Fig. 

6 . Zenith/Brightness: After the bayesian reconstruc¬ 
tion criterion, the remaining air shower background 
events should generally be well described by down¬ 
going reconstructions, although they are still the 
dominant component of the data above the hori¬ 
zon. The air showers are typically dim, however, 
after being attenuated in the overburden, while 
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SUPPL. FIG. 2. The distribution of simulated air shower 
background events and E~'^ signal neutrino events as a func¬ 
tion of reconstructed zenith angle and difference of log likeli¬ 
hoods between reconstructions which do and don’t contain a 
bayesian prior that the events should be reconstructed down¬ 
going according to an approximate, known distribution of 
background event angles. 

the muons produced by neutrinos can be produced 
much closer and so can arrive at the detector with 
much higher energies. By sacrificing all events with 
relatively low energies, as simplistically measured 
by the number of detector modules which receive 
light, it is possible to eliminate virtually all of the 
expected air showers at a given zenith angle, while 
leaving a window for higher energy neutrino events. 
At the same time, very dim events with up-going re¬ 
constructions are also eliminated, as this accounts 
for few neutrinos which are actually valuable to this 
analysis, but also a few remaining air shower events 
with poor reconstructions. This cut is shown in 
Suppl. Fig.§ 

7. Split Event: Despite the separation algorithm ap¬ 
plied during the Data Reduction step, some co¬ 
incident air shower events remain. Therefore, we 


SUPPL. FIG. 3. The distribution of simulated air shower 
background events and E~^ signal neutrino events as a func¬ 
tion of reconstructed zenith angle and number of detector 
channels triggered (NChan). The air showers generally have 
down-going reconstructions and low brightness. 


attempt to remove them by brute force by di¬ 
viding each event into halves, reconstructing both 
halves separately, and testing whether the two sep¬ 
arate reconstructions give a better fit than a sin¬ 
gle reconstruction. Specifically, each event is di¬ 
vided in half both at the mean time of light detec¬ 
tion, and along a plane perpendicular to the sin¬ 
gle particle reconstruction and passing through the 
amplitude-weighted mean position of the light de¬ 
tections. Both division methods are tested, and if 
either yields a better reconstruction the event is re¬ 
jected. The fit is here considered better according 
to the number of photons appearing as unscattered 
with respect to the track hypotheses, so if the split 
reconstructions yield more than twice as many un¬ 
scattered photons as the single reconstruction, we 
conclude that the event is coincident and remove it. 
This final cut eliminates a portion of the remain¬ 
ing air shower background events, but is included 
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mostly for safety against contamination unforeseen 
due to the limited simulation statistics used in de¬ 
veloping the selection. 

Suppl. Tab. 1^ shows the overall results of the data se¬ 
lection process. Background from cosmic ray air showers 
is reduced by a factor of approximately 5.8 x 10^ so that it 
makes up only about 0.1% of the final data, while 23.8% 
of the neutrinos from a hypothetical E~'^ flux which trig¬ 
ger the detector are expected to be retained. 

ENERGY ESTIMATION 

The main power in this analysis to distinguish the sig¬ 
nature of astrophysical neutrinos from the background of 
atmospheric neutrinos is in the different energy spectra 
of these fluxes. However, the neutrino energies cannot 
be measured directly, so in this work the energies of the 
muons produced by the neutrinos are reconstructed, us¬ 
ing the fact that the average energy loss rate of high 
energy muons is proportional to the muon energy. Fur¬ 
thermore, since the muon tracks are not required to be¬ 
gin within the instrumented volume of the detector, the 
energy of the muon on arrival to the detector may be ar¬ 
bitrarily much smaller than its initial energy. These facts 
considerably limit the amount of information which can 
be extracted about the energies of the neutrinos them¬ 
selves. 

In addition, the extracted energy information is also 
limited by the practical capabilities of the muon energy 
reconstruction. The energy reconstruction is selected so 
that the computed proxy has the highest possible res¬ 
olution (among currently available methods), but it is 
not necessarily unbiased and it does not have a one-to- 
one relationship to true muon energies. While an ideal 
reconstruction would have a one-to-one mapping to the 
true physical parameter, fluctuations such as the num¬ 
ber and size of stochastic energy losses, and variations 
of the position of the muon track within the detector 
make this impossible to realize. Bias of the estimator 
can be avoided, or at least largely removed by calibra¬ 
tion, but this serves no purpose in the context of the 
forward-folding maximum likelihood fit, since as long as 
the proxy is related on average to the true parameter by a 
monotonic function, the particular choice of this function 
simply alters the experimental data distribution and the 
simulated template distributions in the same way, and so 
cancels out of the likelihood. In addition, it is not pos¬ 
sible to carry out a calibration procedure fully correctly 
a priori, since the relationship between the true param¬ 
eters and the reconstructed proxy depends on the true 
neutrino energy spectrum. 

While calibration of the energy proxy is not relevant 
for the maximum likelihood fit, the relationship of the 
energy proxy to the relevant physical energies is of gen- 



Muon Energy Proxy [arb. units] 

(a) Distribution of energy proxy values arising from 
different true muon energies. 



Muon Energy Proxy [arb. units] 

(b) Distribution of energy proxy values arising from 
different true neutrino energies. 

SUPPL. FIG. 4. Each column in the figure has been indepen¬ 
dently normalized to form a PDF of possible true parameter 
values. The best fit result with an astrophysical flux 

has been assumed. Fluctuations and missing data at the edges 
of the distributions are due to limited simulation statistics. 


eral interest, and can be explored using the results of 
the fit. Suppl. Fig. shows the results of weighting 
a set of simulated neutrino events to the best-fit spec¬ 
trum produced by the fit and plotting the distributions 
of the muon energy proxy against the true muon energy 
at the point of closest approach to the detector center 
and the true primary neutrino energy. Each bin in the 
energy proxy has been independently normalized, elimi¬ 
nating the influence of the neutrino energy spectrum on 
the distribution of the proxy, making clear the proba¬ 
bility of a proxy value arising from each possible true 
parameter value. The feature which appears at low en¬ 
ergy proxy values (and low true particle energy values) 
is characteristic of the transition to the low muon energy 
region, where energy loss is dominated by ionization and 
varies less strongly with energy, from the high muon en¬ 
ergy region in which stochastic losses dominate and the 
energy loss rate varies more rapidly with energy. For 
energy proxy values larger than ^10^ the most proba- 
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Cut level 

Air shower background 
Rate (Hz) Number in 

full sample 

Conv. atmos 
Rate (Hz) 

neutrinos 
Number in 
full sample 

E ^ signal neutrinos 
Fraction surviving 

Muon Filter 

32.6 

1.86 X 10® 

1.17 X 10~® 

667000 

0.932 

Data Reduction 

10.6 

6.04 X 10“ 

9.13 X 10"® 

520000 

0.897 

Zenith Angle 

3.88 

2.21 X 10“ 

8.08 X 10"^ 

460000 

0.688 

Track Quality 

4.59 X 10"'‘ 

26100 

1.36 X 10~® 

77200 

0.313 

Bayesian Reconstruction 

1.66 X 10"'‘ 

9460 

1.12 X 10"® 

63800 

0.286 

Zenith/Brightness 

2.01 X 10"'" 

115 

5.75 X 10"^ 

32800 

0.238 

Split Event 

5.56 X 10" '' 

31.7 

5.75 X 10~^ 

32800 

0.238 


SUPPL. TABLE I. Rates and fractions of simulated data surviving by type as a function of the level of selection applied. 
Efficiencies are with respect to the detector trigger. Note that the cuts following the ‘Data Reduction’ do not strictly have a 
fully defined ordering; the order shown here is simply chosen for presentation purposes. 



SUPPL. FIG. 5. The flux templates obtained for the best fit 
using a free power law index as functions of primary neutrino 
energy. 


ble true muon energy is somewhat larger (for example, 
^ 2.4 X 10® GeV for a proxy value of 10®). The fac¬ 
tor relating the energy proxy to the most probable asso¬ 
ciated neutrino energy is generally larger, for all proxy 
values, due to the primary neutrino energy being strictly 
greater—sometimes much so—than the muon energy at 
the detector (the most probable neutrino energy leading 
to a proxy value of 10® is ^ 3.3 x 10® GeV). 

Suppl. Fig. shows the energy spectra of the fluxes 
which contribute to the overall best fit. Convolving the 
versions of these distributions and with those of Fig. 4b 
for an E~^ astrophysical neutrino spectrum yield the 
curves in Fig. 2 of the main text. 

















